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Abstract: Gravitational wave production induces a strong constraint on the am- 
plitude of a primordial magnetic field. It has been shown that the nucleosynthesis 
bound for a stochastic gravitational wave background implies that causally generated 
fields cannot have enough power on large scales to provide the seeds necessary for 
the observed magnetic fields in galaxies and clusters, even by the most optimistic 
dynamo amplification. Magnetic fields generated at inflation can have high enough 
amplitude only if their spectrum is very red. Here we show that helicity, which leads 
to an inverse cascade, can mitigate these limits. In particular, we find that helical 
fields generated at the QCD phase transition or at inflation with red spectrum are 
possible seeds for the dynamo. Helical fields generated at the electroweak phase tran- 
sition are instead excluded as seeds at large scales. We also calculate the spectrum 
of gravitational waves generated by helical magnetic fields. 
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1. Introduction 

Magnetic fields are ubiquitous in the Universe. Wherever they can be measured, 
they are found. In stars, in galaxies [|IJ, locally and at high redshift ||, and in 
clusters of galaxies ||. There is also evidence of magnetic fields in super clusters [f|. 
However, the origin of these fields is still unclear. Have they emerged in the late 
Universe from charge separation processes or by ejection from stars and galaxies ||? 
Or have they been amplified from primordial seed fields which may represent a relic 
from the early Universe, from the electroweak (EW) phase transition M or even from 
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inflation |7j? If the second exciting possibility is realized, this means that we can 
learn about processes in the early universe from studying cosmological large scale 
magnetic fields. 

In a previous paper || it has been shown that primordial magnetic fields lead 
to significant production of gravitational waves. If the magnetic field spectrum is 
blue, as it has to be if the production mechanism is causal ||, the nucleosynthesis 
limit for a gravitational wave (GW) background strongly constrains the amplitude 
of magnetic fields on large scales. This strong constraint comes from the fact that 
for causal magnetic fields, the energy density has to behave like 



dlogk 



«r l.i 



with comoving wave number k, on scales which are larger than the correlation scale. 
Hence even the moderate nucleosynthesis limit, since it comes from the smallest 
scales, highest wave numbers, at which the magnetic field is maximal, leads to a 
very strong limit on the field amplitude at large, cosmological scales. The detailed 
results are given in || |T0 . For the derivation of this limit it is assumed that the 



magnetic field spectrum evolves solely via the damping of fields on small scales and 
via flux conservation. On large scales, the magnetic field spectrum scaled to today 
is assumed to remain constant. 

However, if the magnetic field has non-vanishing helicity, the conservation of 
helicity leads to an inverse cascade, i.e. it can move power from small to large scales. 
A derivation of this result can be found in the review \T1\. This can mitigate the 



magnetic field limit which precisely comes from the fact that for causally produced 
magnetic fields there is so little power on large scales. The production of helical 
magnetic fields has been proposed for both, inflation [T2| and the EW phase transition 



where the magnetic field helicity is linked to the baryon number fL3[ . Furthermore, 
the formation of maximally helical magnetic fields at the QCD phase transition has 
been proposed in Ref. |14| . 



The evolution of helical magnetic fields and the inverse cascade have been studied 
in numerical simulations, and simple fits which describe the evolution of the corre- 
lation scale L{t) and of the magnetic field energy density psif) have been derived 
in Refs. |15|, [L6|, [T7[]. Using these results, we want to determine upper bounds on 
the amplitude of helical magnetic fields from the induced GWs. These bounds are 
summarised in Table |1[ 

In Section |^ we present the basic definitions and discuss the evolution of normal 
and helical magnetic fields. Here we make use of the results for the inverse cascade 
discussed in Ref. Jl7|]. In Section ^| we calculate the induced GW spectrum. In 



Section |] we derive the limits on helical magnetic fields on cosmological scales. In 
Section || we conclude. 
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Notation: Throughout this paper we neglect curvature and the cosmological 
constant, which are not relevant for our discussion. The cosmological metric is given 
by 

ds 2 = a 2 (t) (-dt 2 + Sijdtfdx?) , (1.2) 

where t denotes conformal time and the scale factor, normalized to 1 today, is given 
to a good approximation by 



a(t) ~ H t 



90 



9es(t) 



1/6 



1.3) 



where H denotes the present value of the Hubble parameter, g e s(t) is the number 
of effective relativistic degrees of freedom at time t, g = g e s(t ) = 2, and f2 ra d is the 
radiation density parameter today. In the following, the density parameter is defined 
as £lx{t) = px(t)/ p c (t), where p c (t) denotes the critical energy density at time t. 

Spatial vectors are indicated in bold face, 3d spatial indices are lower case Latin 
letters while 4d spacetime indices are lower case Greek letters. 



2. The evolution of helical magnetic fields 
2.1 Basic definitions 

The high conductivity of the cosmic plasma implies that, to lowest order, magnetic 
fields evolve by flux conservation, so that B oc a~~ 2 . We are mainly interested in 
the part of the time dependence of our quantities which is not simply due to the 
expansion of the Universe but to the growth of the magnetic correlation length and 
to the additional decay of the magnetic energy density due to dissipation and to 



the MHD cascade ||T5|, [16|, |17[ . Therefore we eliminate the scaling with redshift by 
expressing all the quantities in terms of comoving ones scaled to today which we 
denote by a tilde. For example the comoving magnetic energy density is given by 


(B 2 (t)) = (B 2 {t))a\t) . (2.1) 

Here (B 2 (t)) depends on time via the evolution of the correlation length and because 
of energy dissipation. Comoving quantities are not multiplied by powers of the scale 
factor a 2 (t) when lowering and rising indices. 



The power spectrum of the magnetic field can be written as ]Tj| 

(S i (k,t)B;(q,t)) = ^5 3 (k - q) [(tfy - kik) 5 (M) + ie ijn k n A(k,t)\ . (2.2) 

The functions S and A denote the parity even and the parity odd parts of the two 
point correlator respectively and k is the comoving wave vector (on the quantities 
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where there is no danger of confusion because they always denote conformal quanti- 
ties like k, S or A, we omit the tilde). Using the above expressions, we compute the 
comoving magnetic energy density ps = (87r) _1 .B 2 , 

f°° dk dp B (k,t) dp B (k,t) _ k 3 S(k,t) 
PB[) ~J k dhgk ' dlogk ~ 2(2vr)3 • [Z - 6) 

Note that the antisymmetric part of the spectrum does not contribute to the energy 
density but its presence, which indicates non-vanishing helicity, influences the time 
dependence of ps and of the magnetic correlation length. In |17| the author derives 
analytically the expressions for the time evolution of the magnetic energy density 
and of the magnetic correlation length. These evolution laws have also been obtained 
in numerical simulations of a magnetic field in a turbulent MHD phase during the 



radiation dominated era 161. On the other hand, the simulations of Ref. \uM indicate 



different exponents for the evolution laws. In the following analysis we adopt the 



analytical picture and the evolution laws derived in ||17|| . It is easy to show that during 
this epoch the MHD equations are identical to those of a non-expanding Universe, 
provided that all physical variables are replaced by comoving variables (see ]n5| Pp. 

Our aim is to compute the GWs generated by an helical magnetic field and to use 
the GW energy density to derive constraints on the magnetic field strength on the 
comoving scale A ~ 0.1 Mpc. For this we express the results of Ref. fl7j in terms of 



the quantities introduced above. Campanelli |17| defines the spectral energy density 
EB(k,t) which is related to Eq. ( |2.3| ) by 1 

dpB{%t) _ k 3 S(k,t) _ 

~dhgir-^(2^- k£B{Kth {2A) 

We adopt the following power spectrum for the magnetic field energy: 

( nu(t) r K < L(t) 

I W) {l+K2)(7+2n)/A for K < j—^ , 

1° fOT K ^T±k- 

Here K = kL(t) /(2tt), L(t) is the time dependent comoving correlation length which 
we infer from |I7| and Zdiss is the comoving dissipation scale (for smaller scales the 
magnetic power spectrum is exponentially suppressed and we thus can set it to 



zero). We derive the dissipation scale ldi SS (t) in Appendix |A.2| . Eq. Q2.5|) is inspired 
by Eq. (19) of ||20|| , which provides a continuous expression for the turbulent velocity 
spectrum, interpolating between the large and small scale behaviours. We extend 
it here to the case of a magnetic field processed by MHD turbulence. The initial 
power spectrum at large scales, kL <C 1 behaves like k 2 S(k) oc k n+2 , n is the spectral 



is using Heavyside-Lorentz units such that pb = (1/2)B 2 while we are using Gaussian 
units with ps = (l/8n)B 2 . This leads to differences of factors of An in the relative expressions. 
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index introduced in Refs. lfL0| , §|. If the initial correlation length is finite, the power 
spectrum given in Eq. (|2.2|) must be analytic and hence n > 2 must be an even 
integer 0. A scale invariant spectrum corresponds to n — —3. We have chosen the 
form of 6B(k,t) such that it scales like k~ 3 ^ 2 in the inertial range, 2ir/L(t) < k < 
2vr//diss- This behavior corresponds to the Iroshnikov-Kraichnan spectrum, which 
can be realised in fully developed MHD turbulence PTf . Here we could have equally 



chosen the Kolmogorov or Goldreich-Sridhar spectral slopes: our main result does 
not depend significantly on this choice. Note that Ref. ]l7j does not account for the 
presence of the inertial range. 



We adopt here the interpolating formula (|2.5|) in order to avoid joining the two 
asymptotic behaviours, at large and small scales, of the magnetic field spectrum 
up to the energy injection scale L(t). This has been done for example in Ref. [§], 
and it leads to an overestimation of the peak amplitude. Note that we extend the 
formula of Ref. [[2IJ to represent also magnetic fields with red spectra, generated 



during inflation. 

Integrating Eq. ( |2.5| ) over k we obtain the total comoving magnetic field energy 
density, 

~ PB (t) = Pc (t )Q B (t) ^ VB (t) r(g^) 2 > ( 2 ' 6 ) 

where we have set the upper limit of integration to infinity. Hence the function 77b(£) 
reflects the time dependence of the total energy density. In the above expression we 
have introduced 

o (A - _ PB(t)a 4 (t) _ PB(t)a\t) 

"M r J = TTT — 7j~\ — TTT— "rad 

Pc[to) Pc[to) Prad(^o) 

\ 1/3 / ,n / \ 1/3 

90 ^ pB{t) n iad = (-*>) n B (t)n I3d . (2.7) 



9es(t)J Prad(i) \9es(t) 

We restrict to the analysis of magnetic fields in a radiation dominated universe. 
Furthermore, adiabatic expansion implies |22| that the entropy S oc g e s(aT) 3 is 
independent of time, so that p ra d(^) = Prad(^o) a (^)~ 4 [5'o/5 , eff(^)] 1 ^ 3 - At time which 
we define as the time at which turbulence is fully developed (as we shall see in the 
following), one has 

ni S n B ( t .) = %M(?i) 1/ ^%M, (2 . 8 ) 

i 'rad V 90 J i 'rad 

where = g c s{t*)- The comoving magnetic energy density parameter f2s(t) depends 
on time via the dissipation of magnetic energy. 

In the following we will often neglect changes in g e s{t). For example, we neglect 
the fact that the evolution of gravitational radiation is modified, even if the universe is 
radiation dominated, due to the fact that while g e s{t) is evolving the scale factor does 
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not expand like a oc t but somewhat faster. For the EW phase transition (g^/go) 1 ^ 3 — 
3.7. In the magnetic field limits this factor enters at most with power 1/2, which 
yields differences of less than a factor of two. In the amplitude of gravitational waves 
it can enter with a higher power and change it by up to an order of magnitude. But 
this is in any case roughly the precision of the results derived in this work. 

2.2 Direct and inverse cascades 

The main difference between non helical field evolution (which only exhibits direct 
cascade and dissipation on small scales) and helical field evolution (which leads to 
inverse cascade) can be expressed in the time evolution of the correlation scale L(t) 
and the comoving magnetic field energy density ps(t), which we cast in the amplitude 

IB® 0- 

We introduce the normalized conformal time r 

r = ^ , (2-9) 

where t* L denotes the initial eddy turn-over time on the scale of energy injection L*, 
and ti n is the time at which the magnetic field is generated. The eddy turnover time 
is defined as t* L ~ L*/(2vl), where vl is the initial eddy turnover speed and the 
factor 1/2 comes from the fact the L* is the eddy diameter. 

Starting from the asymptotic laws given in ||17|| , we rewrite the time evolution 



of the comoving magnetic energy density and of the comoving magnetic correlation 



length in terms of the normalised conformal time (|2.9| ). However, with respect to [17| 
we add a model of the initial evolution in which the magnetic energy density grows 
continuously from zero to the equipartition value, at which MHD is fully developed. 
We do this because later on we evaluate the GW spectrum generated by the magnetic 
source, and the time continuity of the source does affect the resulting GW spectrum 



(see Ref. [^] and section |3]J). Therefore, we assume that the magnetic field energy 
density is zero (continuous) at t — ti n , r = 0; it then reaches its maximal value p* B 
after a 'switching on' time which we set equal to the characteristic eddy turnover 
time, t* L , namely at t* = t in + t* L . We therefore define U = + t* L as the time 
when turbulence is fully developed, and normalise also the energy injection scale as 
U = L(U), L(t in ) = 0. 

• Selective decay (direct cascade): the evolution laws are in this case 

M*) - j r - 2(n+ 3)/(„ + 5) fyrtzu, r>l, (2 ' 10) 
L{t) = L*r 2/(n+5) , (2.11) 

where we may only consider n > — 2 in order to recover the correct behavior 
with respect to time of the above quantities (i.e. decay for the energy and 
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growth for the correlation scale). The above expressions go continuously to 
zero with r — > and they have the asymptotic behavior found in Ref. |17|| for 
t t m > t* L , t 1. Note that Ref. [|I7[] uses a spectral index p which is related 
to n via p = n + 2. For simplicity, the parameter Kjiss of [17] is set to one, 
corresponding to 7 = 3r 2 [(l + p)/2]/(3 + p)/T 2 [p/2]. 



The energy injection scale is determined by the physical process that gener- 
ates the magnetic field and the turbulence. Generically it can be parametrised 
as a fraction of the horizon at the initial time. Therefore, we introduce the 
small parameter e < 1 defined by 

L* = et in , such that t* L ~ = 7^— tin • (2-12) 

A typical value for causally generated turbulence is e ~ 1CT 2 (see for exam- 
ple [[HJ). The necessary condition to have a turbulent cascade is that t* L < t in , 
i.e. vl > e/2. Eddies of the size of the horizon which move at the speed of light 
are the limiting case. Although it grows, the correlation length never becomes 
larger than the horizon. In fact one has 

f(A mi ,*\ / i \ 2 /( n + 5 ) /+ \(n+3)/(n+5) 

^ = [{2v L f e^] V(n+5) (l - *f\ , (2.13) 

which is smaller than one for all times t > t- m . Indeed, for the initial period 
Un < t < U, the term (1 — t ul /t) 2 ^ n+ ^ — > and it dominates the above 
expression, while for t^> t* > t in the asymptotic behaviour is controlled by the 
last term of the equation which keeps the correlation length smaller than the 
Hubble radius. This is shown in Fig. [j]. 

Even though the correlation length is growing, the spectral energy on a given 
co moving scale k is at best constant. On scales which are larger than the 
correlation scale, kL(t) < 1 

dp B (k,t) _ „ U \n„Tu\\n+3 



d\ogk 



cx p B (t)(kL(t)) n+i . (2.14) 



From Eqs. (|2.10| ) and (|2 . 1 1|) it follows that during selective decay L n+3 (t)p B (t) = 
constant, hence dps/d log k does not evolve on large scales. The same behaviour 
is observed in the free decay of the turbulent velocity field, and is related to 
the constancy in time of Loitsyansky invariant (see for example p4" , |25| ) 



Inverse cascade: 

During inverse cascade we assume that the magnetic field energy and correla- 
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Figure 1: Time evolution of the comoving correlation length L(t) as a function of r for 
an inverse cascade (blue, dashed line) and for a direct cascade with the spectral indexes 
n = 2 (red, solid line) and n = —1.8 (green, dotted line). They are compared with the time 
evolution of the comoving Hubble radius in a radiation dominated background, in units of 
the correlation length L* (magenta, dot-dashed line). 



tion length evolve according to fL7 | 



PB(t) = P* B 



for t in < t < U , T < 1 



T 



-2/3 



for t > U 



t > 1 



L{t) = h B L,r 2 '\ 



(2.15) 
(2.16) 



where Hb is the initial fractional helicity: hs = corresponds to a non-helical 
magnetic field that remains non-helical for all its evolution (for which the above 
scaling relations do not apply), while hs = 1 characterizes a maximally helical 
field. The above equations are again valid only for p = n + 2 > 0. 



The original expression of Ref. [117]] for L(t) contains also a logarithmic term, 
and gives asymptotically a slower growth than Eq. ( |2.16| ), as can be seen in 
Fig. |||. In the following we neglect this logarithmic correction. Similarly, ps(t) 
given by Eq. ( |2.15| ) decays more rapidly than the full expression given in [|17f| , 
due to the same logarithmic correction (see Fig. [|). The limits on the magnetic 
field on large scales obtained using Eqs. (|2.16|) and ( |2.15|) are less stringent than 
those one would obtain using the more accurate expression of Jl?j. Neglecting 
the logarithmic corrections is therefore a conservative assumption. 

Although the correlation length grows faster than in the selective decay phase, 
it never becomes larger than the horizon even in this case (see Fig. |l]). During 
the inverse cascade the product L(t)p B (t) = L*p* B is constant. From Eq. 



- 8 - 



10 



Figure 2: Time evolution of the comoving correlation length L(t) as a function of r during 
inverse cascade as given in Ref. Q ( green dashed line) and neglecting the logarithmic terms 
as in Eq. ( [2.16D (red solid line). 



* cq 




10 



Figure 3: Time evolution of the magnetic energy density as a function of r during inverse 
cascade as given in Ref. jl7| (green dashed line) and neglecting the logarithmic terms as 
in Eq. ( 2.15| ) (red solid line). 



one sees that during the inverse cascade the spectral energy density is growing 
like L n+2 (t) at large scales. 



Eqs. (|2.15|) , (|2.16|) apply only after a first phase during which the system behaves 

by selective decay |17] 



as if the magnetic helicity was zero, 



i.e. 



Only when the 



magnetic field (with initial fractional helicity hs) becomes maximally helical, the 
inverse cascade can start, and Eqs. (|2.15|) and ( |2.16|) apply. In order to find the time 
at which this happens, one matches the product L(f)p B {t) (which is time dependent 
during selective decay), to its constant value during the inverse cascade. This defines 
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the transition time th 



i , e ,-(n+5)/(2n+4) 



(2.17) 



For a maximally helical magnetic field th(h,B = 1) = £*■ In general, for a given 
the second stage takes place for times t > th > £* and lasts until the time tfi n at 
which the turbulent phase ends (c.f. next section). Moreover, in the case of zero 
initial helicity, one has pure direct cascade: t h — > oo when h B (this is true only 
if we restrict the value of the spectral index to be n > —2, which we always do in 
the following). 

2.3 The end of the turbulent phase and the dissipation scale 

The turbulent phase ends when the Reynolds number on the scale of energy injection, 
L(t), becomes of order unity [25]. In Appendix |A.1| we calculate the epoch at which 



turbulence ends for the EW and the QCD phase transitions, as well as for inflation 
with T* ~ 10 14 GeV. The most important result from this calculation is that in all 
cases turbulence lasts for many Hubble times and therefore the source is not short 
lived. This finding and its consequences are the subject of For example, for 
a maximally helical field generated at the EW phase transition we find the final 
temperature T fin ~ 21MeV (note that turbulence ends before nucleosynthesis |L6|). 



In Appendix [A.2| we determine the dissipation scale which is the scale l&g a (t) 
below which energy injection no longer leads to turbulence but is simply dissipated. 
This scale determines the time-dependent UV cutoff of our spectra. The evolution of 
both the correlation length L(T) and the dissipation scale ldi ss (T) for the EW phase 
transition is shown in Fig. [|. The dissipation scale grows faster than the correlation 
length and turbulence terminates roughly when the two curves cross. 

In the evaluation of these scales, we often use the approximation T\jT<i ~ £2/^1 5 
which neglects changes in the number of effective relativistic degrees of freedom. 
Moreover in the following we do not distinguish among the temperatures correspond- 
ing to the initial time t 1VL and to the time at which turbulence is fully developed t*, 
since they are separated by less than one Hubble time. Therefore, we generically in- 
dicate with T* the temperature at which the generation mechanism for the magnetic 
field takes place. 

After the end of turbulence the magnetic field simply stays frozen in the fluid 
at scales larger than the dissipation scale Zdiss(2fin) — £(2fi n ). Eventually other 
dissipation processes, due to radiation viscosity, become active 



3. The GW spectrum 

3.1 Generation of GWs from sources 

In this section we calculate the spectrum of the GWs generated by a helical magnetic 
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0.001 0.01 0.1 1 10 100 



T [GeV] 

Figure 4: The evolution of the comoving magnetic correlation length (red solid line) 
and dissipation length (green dashed line), in the inverse cascade phase for the EW phase 
transition. Both quantities have been normalized with respect the initial value of the 
comoving correlation length L*. 



field. We restrict our analysis to the case of maximally helical fields. This calculation 



has also been performed in Ref . |28) , where however stationarity in time of both the 
source and the GWs has been assumed. This invariance under time translation of 
the source seems to us justified only if the time over which the source is active is 
much shorter than a Hubble time. As we have argued in Section |2| and derived 
in Appendix [A.l| , the Reynolds number remains larger than unity for the scales 
of interest over many Hubble times. Therefore we want to re-calculate the GW 
spectrum without the assumption of stationarity. We shall then compare our results 



with Ref. [28 



The parity invariant part the GW spectrum, which is the part which contributes 



to the energy density, is of the form |L0, 

(Ay(M)/i?„(q,t)) = ^A^k^k-q^l^M) 



4 

Here A4iji n (k) is the GW polarization tensor normalized such that M-ijij 
\h\ 2 {k,t) is related to the GW energy spectrum as follows [KJ |T%f : 



(3.1) 



4 and 



Pewit) 



(hijhij) 



2>2nGa 2 {t) 
dp GW (k,t) k 3 \h\ 2 (k,t) 



dk dp GW (k, t) 
k d log k 



dlogk ' 8(2vr) 3 Ga 2 (t) ' 



(3.2) 



Our definition of the metric perturbations differs by a factor 2 with the one 
in JTOJ [Ej and agrees with ds 2 = a 2 {t) \-dt 2 + (<% + h i j)dx i dxP]. The Fourier 
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space expression of the projection tensor onto the transverse traceless component, 



•Myin(k), is given explicitly e.g. in Ref. [[18 



The evolution equation which governs the generation of GWs in an expanding 



universe is simply, see e.g. |22j] 



hij + 2— hij + k 2 hij = \Q-nGa 2 Iiij , (3-3) 

where the source term is the tensor contribution to the anisotropic stress of the energy 
momentum tensor of the source. In our case these come from the magnetic field and 
we relate the tensor anisotropic stress spectrum to the magnetic field spectrum in 
the next section. In terms of the comoving anisotropic stress 11^ = a 4 !!^ the above 
equation becomes 

JL + 2-hu + k 2 ha = 167rG^ . (3.4) 
a a 2 

For GWs generated by a primordial magnetic field we cannot neglect the expansion 
of the Universe, since this source is active over a period which is much longer than 
a Hubble time. In terms of the rescaled variable hy = ahij, Eq. ( |3.4| ) becomes 

hij + (k 2 - hj = 167rG^ . (3.5) 
In a radiation dominated background with a oc t so that a — 0, 

hij + k 2 hij = IStfC^ . (3.6) 

Since initially hij = hij = 0, the solution of the the above differential equation is 
given by the convolution of the source with the retarded Green function Q(k, t\, t 2 ) = 
sin(/c(t 1 - 1 2 )), 

hij(k, x < x &n ) = [ dy ^^v) gin ^ _ y ^ ^ ^ 

& Jx in a {y) 

where we have introduced x = kt\ and y = kt2- 

We assume that the source is active until the final time tg n at which turbulence 
terminates and the anisotropic stress become negligible. This is not completely 
correct, since after this final time we have no longer an inverse cascade, but the 
magnetic field is frozen in and evolves according to flux freezing. However, this 
is relevant only for GW production at scales which are super-horizon at ig n , and 
therefore it does not affect the peak region of the GW spectrum and the value of 



the integrated energy density, which determines our constraints (see Ref. [10] and 



section |3~3| ). 

Once the source has decayed, GWs are freely propagating. This behavior is 
described by the homogeneous solution of Eq. 



hij(k, t > t fin ) = Aij(k) sxn(kt — kt&n) + Bij(k) cos(kt — kt^ n ) . (3.8) 
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The coefficients Ay and Bij are determined by requiring continuity of and hij at 
t — tfin. Matching ( ft.8| ) to the result from Eq. ( |3.7| ) yields 



^i( k ) = -p - y ^ Vj cos(a; fin - y) 

16ttG nij(k,y) , , . 

Sy(k) = — p— y rfy -^j-sin(i fin - y) . (3.9) 

With Eq. Q3. H ), using the above solution for fry, we obtain for t > tfi n 

|A| 2 (M > «»„) = ^(* 2 + H 2 ) «^-> + <B«B«» 

(* 2 + « 2 ) fiftrtjy r- rfB f rf2 C0S(2 _ u Mk,yz) > 



2a 2 \ k* J J x . n J x . n ' " a(y)a(z 

where we have set y = kti and z = kt?,. Furthermore, we have introduced the 
anisotropic stress unequal time power spectrum, 

(fiy(k, *i)n*-(q, h)) = (27r) 3 5 3 (k - q)n B (fc, h,t 2 ) . (3.11) 

To obtain Eq. ( |3.10| ) we have not only performed an ensemble average, but also aver- 
aged over several periods so that (sin 2 (fct)) = (cos 2 (kt)) = 1/2 and (cos(fct) sin(Atf)) = 
0. At times t at which we can observe a GW with wave number k, the latter must 
be largely sub-horizon so that kt 3> 1. We therefore may neglect the second term 
in the pre-factor [k 2 + H 2 ) ~ k 2 {l + l/{kt) 2 ). Rewriting Eq. ([3.10Q as integral over 
time, we find with (|3.2|) 



dpcMk, t) _ 2G _ k3 f'- ^ f'- ^ _ kt2) U B (k t t ) _ 



dlogfc ?ra 4 (t) J t . n J t . n ' x a(*i)a(t 2 ) 
3.2 Magnetic anisotropic stresses 

According to Eq. ( |3.12| ), in order to determine the GWs produced by a cosmic mag- 
netic field, we need to calculate the unequal time correlator of the tensor type mag- 
netic stress, U.B{k, ti, £2), which sources GWs. By statistical isotropy, the tensor type 
magnetic stress 2-point function has the same tensor structure as the one of GWs, 

(27r) 3 r - ~ - ~ 1 

(IIy(k, t)n, n (q, t)) = ^^^ 3 (k - q) [M ijln (k)Il B (k, t, t) + A ijln {k)U A {k, t, t) . 

(3.13) 

In the above expression II a{^, t, t) is the term of odd parity due to the non- vanishing 
helicity of the magnetic field. It it does not contribute to the GW energy density 
but only to their polarization |jl8fl . The odd parity projection tensor is also given 
0. 
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Following fT8]|, we use Wick's theorem to reduce this four point correlator to the 
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convolution of two 2-point correlators. The ansatz (|2.2j ) then gives for the equal time 
correlator 



U B (k, t, *) = M J d 3 g [(1 + 7 2 )(1 + a 2 )S(q, t)S(\k - q|, t) + 4 7 aA(g, t)A(\k - q|, t)] 

. (3-14) 

where we set M\ = 2/ (47r) 5 , a = k • (k — q) and 7 = k • q. In the case of a maximally 
helical magnetic field, the symmetric and antisymmetric parts of the magnetic field 
spectrum are equal on sub-horizon scales, 

\A(k,t)\ — S(k,t) , U>1. (3.15) 

On super-horizon scales helicity is suppressed (see e.g. ||). In order to account for 
this dependence, we introduce the function E(t, q, |k — q|) in the integral fl3.14Q 



v , , 11 1, J 1 ' for gt > 1 and |k - q|t > 1 , 
E(t,q, k-q) = <^ . (3.16) 

U , otherwise , 



and we set 



U B (k, t, t) ~ M J d 3 q [(1 + 7 2 )(1 + a 2 ) + 4 7 «S(t, q, |k - q|)] S(q, t)S(\k - q\,t) . 

(3.17) 

The integral ( p,17|) for the equal time correlator is evaluated numerically and the 
results is approximated by an analytical fit. More details on this are given in Ap- 
pendix |B.1| . Here we simply present the results for two exemplary values of the 
spectral index, a causal spectrum with n = 2 and a red spectrum with n = —1.8: 

n l(Jf , l )4^c) 1+wl ; 3 ; W6) , / „ »=2, (3,8, 

ii B (K, t) ^ ^L\t)j? B (t) 1 l K { $2»/»> ' " = - LS - (3 ' 19) 



Here 



r(l)r(^) 



0.11 for n = 2 , 
0.04 for n = -li 



3.3 The GW spectrum produced by a maximally helical magnetic field 

Let us now consider a magnetic field with maximal initial helicity Hb — 1, which 
immediately (i.e. at = ti n + t2) develops an inverse cascade. To compute the GWs 
produced by this field, we have to make assumptions about the unequal time corre- 
lator of the anisotropic stress. There are different possibilities which are discussed 
in the literature 
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We consider a completely coherent source, namely a source with deterministic 
time evolution for which the unequal time correlator is just the product of the square 
root of the equal time correlators at the different times, 



fl B (k,ti,t 2 



n B (fc,ti,ti)vnB(fc,t2,t : 



(3.20) 



This is not only the simplest approximation, but the results obtained in this case are 
also quite close to the results from a model with exponential decoherence as discussed 
in Ref. [f26fl . Furthermore, for colliding bubbles where numerical simulations exist, the 



totally coherent approximation is in good agreement with the numerical results [23 



This justifies our hope that this approximation captures the main features of the 
resulting spectrum and, especially, that it gives a good estimate for the total GW 
energy density. Note also that this assumption has usually been made in previous 
works, for a magnetic field which is simply redshifting with the expansion of the 
universe [||, [31] . A comparision of different approximations can be found in Refs 

For a completely coherent source we obtain 

2 



rfpcw(M) _ 2G 
dlogk 7ra 4 (t) 



k 



ttin 



dt' cos(H') 



ll B {k,t') 
a(t') 



+ 



'na / / wn B (M0 

dt' sm(kt')— —- 

«(*0 



(3.21) 



In order to compute the above integrals, we substitute approximations (|3.18|) for 
the anisotropic stresses of a magnetic field with a blue spectrum (n = 2), or ( [3.190 
for a red spectrum (n = —1.8). 

We fix the final time at which the source of GWs ceases to be active as the time 
given by the end of turbulence, when Re (L(T^)) — 1. This corresponds to the 
time at which the inertial range (K > 1) is entirely dissipated, when the dissipation 



scale has grown to reach the correlation length, £diss(Tg n ) ~ L(T^'). This condition 
determines the value of the final temperature at which turbulence terminates 
(in principle, the magnetic field is not damped after this temperature, but simply 
stays frozen in the fluid and keeps on generating GWs; however, here for simplicity 
we restrict to GW production during the turbulent phase, an assumption which, as 
previously mentioned, does not affect our result in a relevant way). 

In Appendix |A. 1| we estimate the final temperature for inverse cascade turbulence 
initiated at different times: 



T, 



(i) 



T, 



fin 
(i) 



fin 
(i) 

fin 



21 MeV for the EW phase transition, 
5 MeV for the QCD phase transition, 
1 GeV for inflation. 



(3.22) 
(3.23) 
(3.24) 
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The final time t9^ corresponding to these temperatures is given by |22 



1 En 

«=i(r«)-o.5(«(r«))- 1/6 -^L. (3.25) 

On the other hand, we know that the dissipation scale grows more rapidly than the 
correlation scale. Therefore, when a given wavelength, smaller than the correlation 
scale but initially larger than the dissipation scale, becomes of the order of the 
dissipation scale 2ir/k ~ l diss (T^ (k)), turbulence is dissipated on this scale and 
the GW source has decayed. This defines a second, fc-dependent final temperature 
lS(k)>T^ g^en by 

! diss (T^(k)) = 2n/k = L*/K*. 



Since T^ n (k) > T~J , the final time of integration for the wave number k is £g n (A;) 
(i) 

''fin 



In Appendix A. 2 we derive analytical expressions for t & ^(k), taking into account 



the time evolution of the dissipation length, see Eqs. QA.16] ) and (A.17). The final 
time tfin is given by 

U n (k) = mm[t^,t^(k)j . (3.26) 

Indeed, for scales smaller than L(t^), hence K(t^) > 1, t^ n {k) is equal to t^(k), 
while larger scales are dissipated only at the end of turbulence, t^. 

The ^-dependence of ta n for magnetic fields generated at the EW phase tran- 
sition and at inflation is plotted in Fig. |5|. The final time starts to decrease for 
small wavelengths, namely around K* = K\ ~ 1CT 3 for the EW transition and at 

= K* ~ 1CT 9 for inflation. This value is given by 

^L(Tft n ) ryrif 2/3 T s* ( % V L % ^ ~ ' ' 



1 = ^i(T fin ) = itj-y^ = K{t^ ~ Kl , (3.27) 



so that 

K* ~ I 

2^ L T H 



/v ; ~ ( ^ZkV /3 . (3.28) 



To calculate the integrals ( [3.21 ) we use Eqs. (|3.18|) and ( |3.19|) for the time evo- 
lution of the magnetic anisotropic stress. More details on the explicit form of ( |3.21|) 
are given in Appendix |B.2| . The GW energy spectrum can be written in the form 



dlogk 



3 A/2 (—^] ^Xgw(^), (3.29) 

\90 / "rad 



which holds for magnetic fields with both a blue and a red spectrum. Here g de- 
notes an average number of relativistic degrees of freedom while the source is active. 
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Figure 5: iT* -behavior of ig n for the EW phase transition (green, solid line) and for 
inflation (blue, dotted line), both normalized with respect t™ . 



We neglect this factor in the following (it enters the limits for the magnetic field 
amplitude only as (g / 'go) 1 / 12 ). The integral Z G w(-K*) determines the spectral shape. 
We have calculated it numerically for both cases: a magnetic field generated at the 
EW phase transition with a blue spectrum, with parameters n = 2, e = 0.01 and 
v\ = 0.2, and one generated at inflation with a red spectrum, setting the parameters 
to n = — 1.8 and e = vl = 1. This means that in the inflationary case, the initial 
stirring scale is set equal to the horizon size, and the eddy turnover time is half the 
initial Hubble time, t* L = t in /2 2 . 

For a causal magnetic field spectrum with n = 2, the GW density parameter 
in units of ((l* B ) 2 /^rad is shown in Fig. || Below the peak frequency, located at 
k ~ 2tt /t* L , the numerical result can be approximated by 

e x Kl for < K, < L/(2n t fin ) , 



dtt GW (k,t ) qK A&e 

"OA/; 



\2 



2" 



L B) 



d log k Q 



rad 



e 2 Kl for L,/(27rt fin ) < K* < Lj(2nQ , (3.30) 



k e 3 Kl /2 for L„/(2ttU) < < L,/t* L . 

A part from the peak, the spectrum shows features at wave numbers corresponding 
to the characteristic times of the source: k ~ t^ 1 and k — t^ 1 (t in is too close to to 

2 Note that in the inflationary case, we could as well have chosen the stirring scale to coincide 
with the horizon at any time. One could argue, in fact, that as soon as a scale enters the horizon, 
it causes a stirring of the cosmic fluid, acting as the stirring scale. We have evaluated the GW 
spectrum also setting L(r) oc 'Hr 1 oc t (instead of oc r 2 / 3 ) and consequently pb(j) oc t -1 to 
maintain the inverse cascade. We did not find an appreciable difference among the two resulting 
GW spectra. 
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Figure 6: The GW energy density spectrum from a causal magnetic field n = 2 generated 
at the EW phase transition with maximal initial helicity hs = 1. The spectrum grows 
like k 3 (blue, dotted line) up to the horizon at the end of the turbulent phase, k ~ tl}. 
Then the slope is given by k 2 (green, dashed line) up to the initial horizon k ~ t^ 1 and 
by \fk (red, solid line) up to the peak at k ~ (2ir)/t* L . Above the peak frequency the GW 
spectrum decays rapidly. For the EW phase transition one has //mHz = AK*. 



be distinguishable). More details on the fit, including the values of the parameters 
6j, are given in Appendix |B.3| . Using the above approximation (|3.30|) , we find that 



the integrated energy density parameter of GWs today is approximately given by 



^gw(^o) 



dk dfi GW (A;, t ) 
k d log k 



6 AT, 




2 x 10 



- 8 m 2 



n 



. (3.31) 



rati 



Here we have neglected the decaying part of the GW spectrum, since the main con- 
tribution to the integrated energy density fi GW comes from the part of the spec- 
trum close to the peak k ~ 2ir/t* L . For the numerical value, we have inserted 
e = 0.01 and v\ = 0.2 for a magnetic field generated at the EW phase transi- 
tion with T* = 100 GeV. The GW spectrum from a magnetic field generated at the 
QCD phase transition, T* = 100 MeV, for the same values of e and vl is very similar, 
and in particular it has the same amplitude at the peak. 



Next we consider an inflationary magnetic field with red spectrum n = —1.8. 
In this case we have to use the anisotropic stress given in Eq. (|3.19|) . The resulting 
spectrum is plotted in Fig. ^, and below the peak frequency it can be approximated 




Figure 7: GW energy density spectrum from an inflationary magnetic field with red 
spectrum n = —1.8 and maximal initial helicity Hb = 1. The analytically expected behavior 
is also indicated (c.f. Eq. fl3.32|) ): tx K% n+e (blue, dotted line), oc ^ 2n + 10 )/ 3 (g re en, dashed 
line) and oc ( 2+6n )/ 5 (red, solid line). For inflation one has //Hz ~ 10 7 K*. 



by 



(e 4 ie +6 , forO<^<^, 



cK2 GW (Mo) _, 3A/ -(^) 2 



d log k Q 



rati 



e6 /C (2+6 " ,/5 , for (<1J1) 1/7 <A'.<# 



L 



Here, the if*— dependence is written in terms of the general spectral index n, and it 
is valid for any n < —3/2. On the other hand, the values of the matching constants 
€i are derived in Appendix p.3| under the assumption n = —1.8. 

Neglecting the decaying part of the GW spectrum, and using the slope j^~( 2+6n )/ 5 
up to the peak k ~ 2ir/t* L , we find the total GW energy density parameter: 

~ * 2 / ~ \ ~( 2+6ra V 5 

^Gw(to) * 3 Af 2 ^- ( _Z|_ ) e6 [ ^ ) n < _3 (3>33) 



(,o*^ 2 

5.2 



rati 



Note that this approximation causes an overestimation of the total GW energy den- 
sity of about three orders of magnitude. However, this does not affect the bounds on 
the magnetic field amplitude significantly: it translates into a bound that is stronger 
by about the 20% (see section |j). For the numerical value in the last equality of 
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Figure 8: The GW amplitude h(f) as a function of frequency from a maximally helical 
magnetic field with £l* B ~ 0.1. Blue, dotted: causal generation at the EW phase transition 
with n = 2; red, solid: generation during inflation with n = —1.8. 



([3.33 ), we have inserted the value n = —1.8 and e = v l — 1. 



Finally, to make contact with future observations, we express the GW spectra 
in terms of the GW amplitude as function of the frequency /. For this we use |32] 



h(f) = 1.26 x 10- ls h y/n GW (f) (yj , (3.34) 



where the frequency is / = k/(2ir) and 

df2 G w(Mo) 



Ogw(/) 



dlog k 



k=2wf 



The behavior of h(f) for a causal and for an inflationary produced magnetic field, 
choosing a maximal magnetic field amplitude of Q* B ~ 0.1, is plotted in Fig. || 
For another magnetic field density parameter the resulting amplitude h(f) is simply 
rescaled by the factor Q* B /0.1. 

The result for a magnetic field with a blue spectrum, shown in Fig. 13, agrees 



qualitatively with the one obtained in Ref. []28j . The position of the peak is the same, 



although the amplitude is higher in Ref. [[28| by nearly one order of magnitude. This 



may come from the different assumptions for the unequal time correlator: in [28| the 
source is assumed to be stationary with exponential decorrelation, while we choose 
a completely coherent source. Moreover, in the magnetic energy spectrum, we have 
modelled the transition from the k n+2 — behaviour at low wave number to the inertial 
range, while Ref. extends the Kolmogorov spectrum down to k — 2tt/L # . We have 
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experienced that this leads to a significant overestimation of the peak amplitude [M] 



Furthermore, the low frequency tail of the spectrum in [^8 grows as y/f, whereas in 



our case it becomes constant for / > l/tfi n : this is due to the fact that our source is 
long lasting, while theirs lasts for less than one Hubble time (see [26|). 



For Q* B ~ 0.1, the EW result is somewhat below the sensitivity range of LISA 
[33|j . On the other hand, the inflationary result is much below any proposed exper- 



iment: even though the energy density is higher, much of it is at high frequency, 
resulting in a very low gravitational wave amplitude (see Fig. |8[). 

3.4 The peak position 



In Ref |2l| it has been argued that approximating the source of GWs by a discon- 
tinuous function might influence the decay law and even the peak position of the 
resulting GW spectrum. Therefore, as already discussed in section [2.2| , in this work 
we model the 'switching on' process and avoid a discontinuous source function. Here 
we compare our result with what we would have obtained assuming a discontinuous 
source. 

Fig. |9| shows the results for a continuous source and for a discontinuous one where 



the inverse cascade starts instantaneously at t = £*. The predictions of Ref. 
are confirmed: in the discontinuous case the peak is no longer at t* L but rather at 
~ t* L /(2vL), leading to an over estimate of the resulting GW output. Therefore, 
it is important to take into account continuity. 

In order to further clarify this issue, we vary the initial speed vl, which relates 
the initial correlation length L* and the characteristic turnover time t* L ~ L*/(2vl)- 
We fix it to vl = 1CT 2 so that and t* L are clearly separated: t* L ~ 10 2 !/*. As can 
be seen in Fig. || in the discontinuous case the peak position is independent of the 
velocity vl, while in the continuous one it is located at k = 2%/t* L . 

Summarising, having assumed that the magnetic field processed by MHD tur- 
bulence needs a characteristic time of order t* L to 'form' induces a peak in the GW 
spectrum at a wave number corresponding to this characteristic time. If the field 
builds up much faster, almost instantaneously, the peak can move to k = 2tt/L*. 



4. Limits 

The main aim of this paper is to derive constraints on the amplitude of a primordial 
magnetic field scaled to today, smoothed over an interesting characteristic scale. We 
choose the comoving scale A ~ 0.1 Mpc. On one hand, this is approximately the 
largest correlation scale of cluster magnetic fields today, on the other hand it is only a 
little larger than the smallest scale which survived dissipation prior to recombination 



2%. Note however that, accounting for the full evolution of the magnetic field, 



Ref. |TB[ found a smaller dissipation scale at recombination, of about 1 kpc. 



- 21 - 




Figure 9: The GW spectrum, normalised to the magnetic energy density, for a causal 
magnetic field n = 2 with maximal initial helicity hs = 1. Continuous source: green, 
dashed line with v\ = 0.2 and magenta, solid line with v\ = 10~ 4 . Discontinuous source: 
orange, dashed line with v\ = 0.2 and blue, solid line with v\ = 10 -4 . In the continuous 
case the spectrum peaks at k ~ 2ix jt* L while in the discontinuous one it peaks at k ~ 27r/L*. 

To constrain the magnetic field amplitude, we use the GW energy density gener- 
ated by the magnetic field and we apply the nucleosynthesis bound to these GWs, the 
same strategy followed in Ref. ||. We define the comoving magnetic field smoothed 
on the comoving scale A by 

4(A) = yj d 3 xB^)exp , (4.1) 

where V\ is the normalization volume given by 

V, = J d 3 yexp f-t) = (v^A) 3 , (4.2) 

and B is the magnetic field scaled to today, -B(x, t) = B{x) / a 2 {t). A short calculation 
gives the smoothed amplitude 

?2 1 f tW__ , ^ ( y 2 



X= V 2 J d x{ PW ■ B(x + y)> exp , (4.3) 

whit V2 = (V2rr\) 3 . Translating the above expression to Fourier space we obtain 

I f°° I k 2 \ 2 \ f°° ( k 2 \ 2 
5^ = —/ dkk 2 S(k,t)expi — =8tt/ dke B (k,t)expi — 

_ ' (4.4) 

We relate this amplitude to the comoving magnetic density parameter Q* B using 
Eqs. (|2.5|), (|2.6|) for SB{k,t) which are valid for both direct and inverse cascade. 
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Performing the above integral we find 



ffi(t) = 8np c (t )Q B (t) r I 1 U 




(4.5) 



where U denotes the confluent hyper-geometric function, see e.g. p4] . 

We assume that the initial helicity of the cosmic magnetic field responsible for 
GW production is maximal, hs = 1. Therefore, during the inverse cascade phase, the 
magnetic correlation length evolves as given in Eq. (|2.16 ), and the product L(t)fle(t) 
is constant in time. The inverse cascade goes on until the temperature T^ n at which 
MHD turbulence terminates. In the maximally helical case, this can vary from 
a temperature of a few MeV to approximately 1 GeV depending on the epoch of 
generation of the field (see Eqs. ( |3.22| )-( |3.24| )). Up to this final time, the correlation 
length has grown substantially, but one readily confirms that it remains several orders 
of magnitude smaller than our scale of interest, A = 0.1 Mpc ~ 10 13 sec, 



L(t 



fin/ 



L*. 



2/3 



< A 



(4.6) 



Therefore, we can expand the function U(a,b,z) = z a [l + 0(z With this we 
find 



n+3 



(V27r)«+3r(i)' 



A 



(4.7) 



Taking into account that during inverse cascade fis(t)L(t) = constant and L(t) 
L^t 2 / 3 , we obtain at t* <^ t < £g n 



Bl(t) 



Bl{U)r* n+ W 3 , with 

u^F^) Pcio) B It 



n+3 



(4.8) 
(4.9) 



For t > tfi n the primordial fluid enters in the viscous phase and the magnetic 
field energy density is dissipated by radiation viscosity |2"7] . During this phase, the 
evolution on large scales is similar to MHD selective decay, i.e. such that the large 
scale part of the power spectrum remains constant: f2s(i).L n+3 (t) = constant. From 
the general evolution of B\ given in Eq. (|4.7p , we see that it is justified to assume 
that, on the scale A, after t nn the magnetic field energy density evolves only by 
redshifting. Therefore, the comoving quantity B\(t) remains constant: 



Bl 



B&t > tfin) 



BKQ [ l f L 



2(n+2)/3 



(4.10) 
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Like every contribution to radiation energy density prior to nucleosynthesis, f^cw 
is constrained by the nucleosynthesis bound |35| 



^GW < ftlim - 0.1 ft Iad . (4.11) 

Via Eqs. (|3.31|) , ( |3.33| ) this yields a constraint on the magnetic field energy density 
parameter fl* B , in terms of Z/*/ti n = e and L*/t* L ~ 2v L : 

& B < (ell g )1/12 ^ad blue case n = 2, EW , (4.12) 
tt* B < 0.144 2+6n)/1 °^rad red case n = -1.8, inflation. (4.13) 



With the help of Eqs. ( t4.8|) -( |O0l) we translate this into a constraint on B\. Using 



Prad(io) = Pc(t ) firad ^ 2 x l(T 51 GeV 4 ~ 0.4 x 1(T 12 (Gauss) 2 , we find for a blue 
magnetic field generated at the EW phase transition 



B\ , 0.3 /M ( " +3,/ V^"' +2,/3 



< 



(4.14) 



1/xGauss ~ {e n v 2 L y/ 24 I A J \t* L 

for n = 2, generated at the EW phase transition 

and for a red magnetic field generated during inflation 

, v (n+3)/2 N ( n +2)/3 
5^ ( 2+6n)/20 /taA ( } 

1/iGauss ~ L ^Ay V *l / 

for n = —1.8, generated at inflation. 

In these equations the pre-f actors are calculated using the spectral indexes and the 
initial and final temperatures corresponding to the generation times. The dependence 
onj/i and e is kept explicit for completeness. 

It is interesting to see that the inverse cascade simply relaxes the limit by the 
factor Tg^ +2 ^ 3 , absent in the non-helical case for which B\(t) = B\. This factor 
tends to 1 for n — > —2, the limiting value for which the inverse cascade relations 
( [2. 15 ) and ( 2.16|) apply. For causal generation with n = 2 the limit for a magnetic 



field is substantially reduced, while for a red magnetic field spectrum with n = —1.8 

2 /3 

the reduction is only by r fin . 

Let us apply our findings to the two generation mechanisms considered above 
(EW phase transition and inflation), to which we add also the interesting case of the 



QCD phase transition []rjj. As mentioned in Sec. |3.3| , we have evaluated the spectrum 
also in this case, finding a very similar amplitude to the EW phase transition case, 
for the same values of vl and e. Therefore, we are confident that we can trivially 



extend the above Eq. ( [1.14j ) also to this case. 



For a maximally helical field we find the following limits on B\. 
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• If the field is generated at the EW phase transition at 100 GeV, t [n ~ 7.8 x 10 4 
sec, assuming a causal spectrum with n — 2, taking the values e = 0.01, 
v\ = 0.2, T fin ~ 21 MeV, setting c/ eff (T fin ) = 43/4 and using A = 0.1 Mpc 
~ 10 13 sec, we obtain the constraint 

^o.iMpc < 8 x 10~ 24 Gauss, EW phase transition, % = lOOGeV (4.16) 

• If the field is generated at the QCD phase transition at 100 MeV, t- in ~ 1.1 x 10 8 
sec, with the same parameters as before but Tg n ~ 5 MeV, the constraint 
becomes 

5 .i Mpc < 2 x 10~ 19 Gauss, QCD phase transition, T* = lOOMeV (4.17) 

• If the field is generated during inflation at T* ~ 10 14 GeV, t in ~ 7 x 10~ 8 sec, 
with an acausal red spectrum n = —1.8, choosing e = vl = 1, ~ 1 GeV 
and (?eff(^fin) = 287/4, we find the constraint 

5 .iM P c < 2 x 10~ 18 Gauss, inflation, n = -1.8 (4.18) 

• If the field is generated during inflation but with a blue, acausal spectrum 
n = 0, with the same values of the parameters as before we find the stronger 
constraint 

^o.i Mpc <4x 10" 28 Gauss, inflation, n = . (4.19) 

The above limits are summarised in Fig. [10] as a function of n, and in Table |I[ 
These upper bounds on the amplitude of the primordial magnetic field are less strin- 
gent than the ones obtained from a direct cascade by the factor t^ +2 ^ 3 . Moreover, 
they strongly depend on the choice for the smoothing scale A, in particular for blue 
spectra. The scaling with A is in fact given by (c.f. Eq. ( |4.7| ) and 0) 

(n+3)/2 

B Xl = Bx 2 \'-~ ± \ . (4.20) 




Therefore, for a smaller smoothing scale of e.g. A = 1 kpc |16|| , the above bounds 
are relaxed by a factor of 10 5 in the EW and QCD generation cases, and by a factor 
of 10 3 in the inflationary case with flat spectrum n = 0. For red spectra the bound 
does not change much, e.g. in the inflationary case with n = —1.8 it is relaxed only 
by a factor of about 16. 

More stringent bounds on the amplitude of a causally produced magnetic field 
can be obtained by imposing simply that the energy density of the magnetic field 
cannot overcome 10% of the total energy density in radiation at generation time: 

^<o.in rad . 
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Figure 10: Upper bounds on the comoving amplitude of a primordial magnetic field from 
GW production as a function of n, for A = 0.1 Mpc, for a field generated at inflation (blue, 
dotted), and at a phase transition (EW: red, dashed, QCD: green, solid). In these last two 
cases, the generation is causal, consequently only n > 2 is allowed. 



Comparing this last inequality with Eq. ( |4.12| ), one sees that in the latter the factor 



multiplying f2 ra< j is 5 orders of magnitude larger than 0.1 (with the usual values 
for the parameters). This shows that, in the causal case, the conversion of magnetic 
energy density into GW energy density, although quite efficient, is not at all complete. 
For example, using Q* B < 0.1f2 ra d, the bound on an helical magnetic field becomes 
-Bo.iMpc ^ 5 x 10 -26 Gauss for the EW phase transition. In the inflationary case, on 
the other hand, the bounds are not modified, since the conversion into GW is much 
more efficient (c.f. the pre-factor in Eq. ( |4 . 1 3|) ) . 

However, accounting for GW production seems to us more model independent. 
Once GWs are generated they do not interact with the cosmic fluid and simply 
redshift with the evolution of the universe. We are therefore sure that any GW 
energy density sourced before nucleosynthesis is still present at that time and must 
respect the nucleosynthesis bound. On the other hand, magnetic energy density can 
be dissipated or converted into other forms of energy during the evolution of the 
universe. We could therefore invoke the extreme scenario in which a magnetic field is 
formed in the very early universe with amplitude higher than the presumed radiation 
energy density at that time. The magnetic energy can subsequently be transformed 
into other forms of energy in such a way that it satisfies the nucleosynthesis bound 
at nucleosynthesis. The main motivation to consider GW production is to obtain a 
bound which is safe from these exotic, but in principle possible, scenarios. 
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5. Conclusions 



In this paper we have derived new upper bounds on the amplitude of primordial 
magnetic fields. We have considered helical magnetic field power spectra, which 
evolve via inverse cascade transferring power from small to large scales. 

For the non-helical case, upper bounds on the magnetic field amplitude on the 
cosmologically relevant scale A 0.1 Mpc have been derived in previous analyses 
[§, These bounds apply to magnetic fields generated before nucleosynthesis, for 



which the correlation scale at the moment of generation is <C A. If the magnetic 
field spectrum is blue n > —3/2, the peak of the energy density per logarithmic scale 
sits at L*, then the amplitude at the scale A is very constrained: 



On the other hand, if an inverse cascade is active, the power at L* is moved to 
the larger correlation scale L(r), following the evolution law L{r) = L^r 2 ^ 3 . At the 
end of the inverse cascade process, we have seen that the magnetic field on the scale 
A is finally 



(n+3)/2,_>+2)/3 



fin 



(5.1) 



(naively one might expect a scaling like t^ +3 ^ 3 , but some of the initial amplitude is 
lost during the inverse cascade process, so that ClB(i)L(i) = constant. This reduces 
the growth of B\ by a factor r fi ^ ). 

The strong limits from magnetic fields which obey a direct cascade are therefore 
mitigated in the helical case by a factor 



(n+2)/3 



fin 



2^l T* / 



1/6" 



(n+2)/3 



e Tfin VS'fm, 

5.3 x 10 7 for the EW transition, 

n = 2 , % — lOOGeV 
2 x 10 4 for the QCD transition, 

n = 2 , T* = lOOMeV 
9.3 for inflation, 

n = -1.8 , T # = 10 14 GeV 
4.8 x 10 9 for inflation, 

n = , % = 10 14 GeV . 



(5.2) 



In this paper we have only considered spectral indexes n > — 2, since Ref. [17] does 
not analyse smaller spectral indexes, and numerical simulations have always chosen 
either n = 2 fll5]| , or n = (T^|. Consequently, we do not know whether for a red 
spectrum with n < —2 the inverse cascade is still active, or whether the limits are 
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those of the direct cascade. For spectral indexes close ton = —2 the above mitigating 
factor is small, even if the magnetic field is generated in the very early universe. On 
the other hand, for blue magnetic fields generated at the QCD phase transition, 
the inverse cascade is not very efficient, since turbulence anyway stops after e + e~ 
annihilation: therefore T^ n is not very large. 

In Table [TJ we summarise our results. We give the upper bounds on the magnetic 
field amplitude obtained both accounting for GW production and imposing that the 
magnetic energy does not overcome 10% of the radiation energy, for helical and non- 
helical magnetic fields smoothed on scales of 0.1 Mpc and 1 kpc. The upper bounds 



for a non-helical magnetic field are in agreement with those given in [p|, |10|. 

We have found that only red magnetic field spectra from inflation or helical fields 
from the QCD phase transition can have the amplitude of -£>o.im P c ~ 10~ 22 Gauss 
which is necessary for amplification by a dynamo mechanism up to the observed 



/iGauss field |TT]. Especially, the well motivated helical fields from the EW phase 
transition are still too constrained, even after the inverse cascade. 

This leads us to the conclusion that the observed magnetic fields in galaxies and 
clusters have either not been seeded by primordial fields, or these primordial fields 
have been produced during inflation and have a red spectrum n <, — 1.8, or they 
have been produced during the QCD phase transition. In this latter case it is crucial 
that these QCD fields be helical because the boost by the factor of about 20000 is 
absolutely needed, while for red (—2 < n < —3/2) inflationary fields the inverse 
cascade is not relevant. 

To evade this conclusion one can argue that magnetic fields coherent on a smaller 
scale, of about 1 kpc, and with the required amplitude of 10~ 22 Gauss are sufficient 
to give rise to the fields observed today in galaxies and clusters ||16| . If this is so, then 



the bounds derived here are relaxed in such a way that also helical fields from the 
EW phase transition can have a sufficiently high amplitude (however, for non-helical 
fields this is still not enough - see Table p. 
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Appendix 

A. The end of turbulence and the dissipation scale 
A.l The end of turbulence 

The turbulent phase ends when the Reynolds number on the scale of energy injection 
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GW limits 




helical 


non-helical 




A = O.lMpc 


A = Ikpc 


A = O.lMpc 


A = Ikpc 


EW 
QCD 

Infl. n = -1.8 
Infl. n = 


8 x 10~ 24 G 
2 x 10 19 G 
2 x 10 18 G 

4 x 1(T 28 G 


8 x 10 19 G 

2 x 10 14 G 

3 x 10 17 G 

4 x 10~ 25 G 


10- 31 G 
10- 23 G 
2 x 10 19 G 

8 x 10- 38 G 


10- 26 G 
10 18 G 
3 x 10 18 G 

8 x 10~ 35 G 


limits from fl* B < 0.1 


EW 
QCD 

Infl. n = -1.8 
Infl. n = 


5 x 1(T 26 G 
10 21 G 
2 x 10 18 G 

4 x 10- 28 G 


5 x 10 21 G 
10 16 G 

3 x 10 17 G 

4 x 1(T 25 G 


6 x 10- 34 G 
6 x 10~ 26 G 
2 x 10 19 G 

8 x 10- 38 G 


6 x 10~ 29 G 
6 x 10 21 G 
3 x 10 18 G 

8 x 10~ 35 G 



Table 1: This table summarises the upper bounds for the magnetic held amplitude 
averaged over the scales A = 0.1 Mpc and A = 1 kpc, for the different generation epochs 
discussed in the paper. Here we present the limits for maximally helical as well as non- 
helical helds. In the four top rows we give the limits from the production of GWs while 
in the lower part of the table we present the limits coming from the requirement that the 
magnetic field contribution be always subdominant: more precisely we require £l* B < 0.1. 
The values which may be sufficient for dynamo amplification are given in boldface. 



becomes of order unity [B^j 

v L L 



Re(L(T fin ),T f 



fin J 



V 



~1. (A.l) 



Here L denotes the physical energy injection scale, v is the kinetic viscosity, and vl 
is the eddy velocity on the scale L. We assume that in the MHD cascade kinetic 
and magnetic energy have the same time evolution (equipartition). Substituting 
L(t) = L(t)a(t) = L(t) (T /T(t)) (g / g eS (T)) 1 / 3 7 we find (<7 eff (T*) = g m ), 

Re(L(T),T) = R^(-^) ^r a , = Re (L*, T*) , (A.2) 

where the power a represents the evolution of the product L(t)vL(t) and we use 
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v L (t) oc y/f%. With Eqs. (|0| pTD and flTIS| , gig ) we obtain 

^ _ f — (n + l)/(n + 5) for non- helical fields (normal cascade), ^ ^ 

\ 1 /3 for helical fields (inverse cascade) . 

The kinetic viscosity is approximately given by the mean free path of the particle 
with the weakest interaction [30], v ~ £ m f p /5. During the early radiation dominated 



phase, neutrinos determine the viscosity until they decouple at T ~ 1 MeV, when 
photon viscosity sets in. For T > 1 MeV the mean free path of the neutrinos is given 
by 

il^j^E, lMeV<T<100GeV, (A.4) 

where Gf = (293GeV)~ 2 is the Fermi coupling constant. Below 1 MeV we have 
to consider the photon mean free path which can be approximated by Thomson 
scattering 



*M ^ 1 3m 2 e m p /V s 

' mfp a T n e 87re 4 tt bPc \ T 



0.3 eV < T < 1 MeV , (A.5) 



where <jt is the Thomson cross section, n e is the electron density, Q b is the baryon 
density parameter and we neglect the short period of time during which electrons are 
still relativistic, after neutrino decoupling but before electron positron annihilation. 
The lower limit in ( |A.5| ) comes from recombination, when the electron density drops 
sharply and photons decouple. The situation also changes at very high temperature, 
when the EW symmetry is restored, T > lOOGeV. Then the coupling constant is 
nearly independent of temperature, and the relativistic mean free path is of the 



order of the inverse temperature |36 



22 

^mfp — j T > 100 GeV . (A.6) 

We first calculate the Reynolds number at the time of generation and at the 
initial correlation scale of the magnetic field, in order to confirm that a turbulent 
MHD phase is indeed present. With (|A.2|) we obtain 

, rp / \ 1/3 

'in J / 90 



R*=v L e-^i^ , (A.7) 

where we have identified the temperatures T in ~ T*. Furthermore, we use v = £ m f p /5 
and 

1 -i/6 rapi 



^in = r— ~ O' 5 ^* 







R.~*aid 2 mpi ( A. 

y/ 9* -^*^mfp(^*) 
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We consider the two situations, 1 MeV < T* < 100 GeV and T* > 100 GeV, where 
the expressions (|A.4|) and ( |A.6|) for the mean free path give 



~ 2stJ mpi 7*Gj. , 1 MeV < T, < 100 GeV (A.9) 



R *~ i^H^' T„ > 100 GeV . (A.10) 

Let us start by considering the generation of turbulence during the EW phase 
transition at 100 GeV. Setting % = 100 GeV, e = 0.01 and v\ = 0.2, with g* ~ 100 
we obtain ~ 10 13 . The corresponding parameters for the QCD transition at 
T* = 100 MeV and ~ 10 yield R* ~ 10 4 . Both Reynolds numbers are much larger 
than one so that we can be certain that a first order phase transition will induce 
turbulence. To determine when turbulence terminates we use Eqs. (|A.2j ) and (|A.4j) , 
setting 

1 - Re^),^) - ^ (^-y ^ (f )' , R , ffi . 

(A.ll) 

For the last equal sign we have used t* L = L*/{2vl) = eti n /(2vi), t ~ t/t* L and we 
have approximated tfin/^in — T^/T^ n . This corresponds to neglecting changes in the 
number of effective relativistic degrees of freedom. 

For direct cascade with n = 2, hence a = —3/7, we obtain T^ n ~ 200 MeV for 
the EW phase transition and Tg n ~ 20 MeV for the QCD phase transition. 

In the helical case with inverse cascade a = 1/3, turbulence is maintained longer 
and we find Tfi n ~ 21 MeV for the EW phase transition and Tfi n ~ 5 MeV for the 
QCD phase transition. 

For generation of magnetic fields and turbulence at the end of inflation T* = 
10 14 GeV, no causality restriction holds and we choose e — Vl — 1. As long as 
v{T) (x £ m f p ~ 22/T, the Reynolds number at the correlation length L(T) evolves like 
(t/t in ) a ~ (T*/T) Q according to Eq. Q (where we neglect changes in the effective 
number of relativistic degrees of freedom). At T* ~ 10 14 GeV, with g* ~ 200 we have 
R* ~ O.Olmpi/T* ~ 10 3 . As time evolves, the Reynolds number at L(T) decays only 
in the non-helical case if n > — 1, so that a = — (n + l)/{n + 5) < 0. In the helical 
case and for n < — 1 the Reynolds number Re(L(T),T) grows as the temperature 
drops. Once T = T ew = 100 GeV is reached, the viscosity v starts decaying rapidly, 
like T~ 5 , and the Reynolds number then decreases. 

We consider three cases 
i) Direct cascade with n = 0, a = — 1/5 



Re (L(T ew ),T ew ) = R*(TjT ew ) a ~ 6 



so T fin ~ 100 GeV. 

ii) Direct cascade with n = —3/2, a = 1/7 
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In this case 

Re {L{T ew ),T ew ) = R*(T*/T ew ) a ~ 9 x 10 4 , 

and 

T fin = T ew Re (L(T ew ),T ew )^ ~ 5 GeV . 

iii) Inverse cascade, cc = 1/3 

In this case 

Re(L(T ew ),T ew ) = R*(TjT ew ) a ~ 2 x 10 7 , 

and 

T fin = T ew Re (L(T ew ), T ew )^ ~ 1 GeV . 

We draw the important conclusion that in all cases the MHD turbulent phase 
always lasts for many Hubble times before the total kinetic energy is dissipated [^(j] . 

A. 2 The dissipation scale 

In the previous subsection we have considered the energy injection scale L(T) and 
determined first that turbulence is present on this scale, and second when turbulence 
ends (i.e. when the entire Kolmogorov range is dissipated). Now we want to know, for 
a given fixed time t (or temperature T), what is the scale below which kinetic energy 
is dissipated. This defines the comoving dissipation scale / d i SS (T). The function 
l<Hss{T) can be found considering that, on scales smaller than this scale, viscosity 
dominates, therefore there is no turbulence. Thus, /diss corresponds to the physical 
scale at which the Reynolds number is equal to 1, 

Re(/ diss ,T) = ^~l. (A. 12) 

v 

Here v is the kinetic viscosity as in the previous section, and v\ is the eddy velocity 
at the dissipation scale. We determine / d i SS (T) only for the helical case which is our 
main interest in this paper (see p6| for the non-helical one). In the inertial range 



the turbulent eddy velocity obeys a Kolmogorov spectrum so that JTU 



Now we use Re(L(T),T) = vl(T)L(T) / V(T) so that, from Eq. ( |A.2j ) neglecting 
changes in the number of relativistic degrees of freedom, we find 

Re ( id ,„ T) = Re (L(T),T) K ^ 

Hence Re (/diss, T) = 1 yields 

/diss(T)^/\ s r 5 / 12 (|^p) 3/4 , (A.15) 
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where we define Z^ iss = L*/i2* . For the last equal sign we have used the behavior 
of the correlations scale with r = (t — ti n )/t* L as r 2 / 3 , (|2.16|) , and a = 1/3 for the 
helical case. 

The evolution of both the correlation length L(T) and the dissipation scale 
l<Hss(T) are compared in Fig. |]for the EW phase transition. Turbulence stops roughly 



(i) 

fin • 



when the two curves cross. We call this time 

Finally, we have to take into account that for a given comoving scale, I = 2tt jh 
the Reynolds number can become unity long before the end of turbulence. The time 

~ (2) 

at which turbulence on the scale I is dissipated is denoted t R ^(k) and it is defined by 



1 = 2ir/k 



'■diss 



Let us first consider Tl < T P , 



100 GeV, so that for all times of interest the 



kinetic viscosity behaves as v oc T 5 . This leads to 



2tt 



*S(*) 



^diss 



1 u 



2vj 



5/12' 



12/41 



(A.16) 



(2) 

Here again, we neglect a possible difference in g e ^ between t* and tg n and we set 



T, 



(2) _ + (2) 



fin 



&n 2v L /(et h 



We also use R, = (L*/l* diss ) 4 / 3 . 



The situation is somewhat more complicated for generation temperatures T* > 
T ew : until T ew the kinematic viscosity decays roughly like v oc 1/T. We therefore 
have to distinguish between scales which are damped at temperatures above T ew and 
those which are damped below. Since we are in this situation only for the inflationary 
case, we set e = vl = 1 for this case. We then obtain 



«(*) 



2 V * 1% 

diss 



12/5 



for such that 4n (^) < *ew , > K* e 



12/41 



2-5/41 for K < r: _ 



(A. 17) 



Here K* is the value of if* for which turbulence terminates at TL 



K* ~ 

ew 



T 



2T* 



5/12 



/* 

Miss 



(A.18) 



When applying these formulas for the calculation of GWs, we must choose the 
true final time given by 

t Rn = mm[4lt^(k)] . (A. 19) 

For turbulence from the EW phase transition and from inflation this function is 
plotted in Fig |] . At turbulence is dissipated on all scales I < L(t^) which is 
the scale of the largest eddies: therefore, the entire Kolmogorov range is dissipated. 
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B. The equal time correlator and other integrals 



B.l U B (k,t) 

According to Eq. ( p.lTf ) we can write the equal time correlator as the following 
integral 



n B (M) 

h(k,t) 



-±L\t)pl{t) [h(k,t)+AI 2 (k,t)} = -±L\t)~p 2 B {t)I{k,t) , (B.l) 



2tt 



2tt~ 



L(t)/l diss (t) 



dQ- 



Q' 



X 



1 + Q2)(2n+7)/4 

rf 7 (l + 7 2 )" 



— 4— 



£2 

2AQ; 2 .. 
dis 



I + x 2)(2n+7)/4 



I 2 (M) 



[2^ 2 + (i + 7 2 )g 2 

£(t)/'dias(t) 

L(i)/(2ni) 



Q n + 2 



(B.2) 



-i 



2Q IK 



1 + g2)(2n+7)/4 

c?7 



x 



f-' 



2JfQ(27rt)' 



n-l 



2KQK. 

diss 



( 1 + x 2)(2n+7)/4 



7 (K - Q 7 ) 



(B.3) 



where we have used a = (k — 7 g) /\Zk 2 + q 2 — 2^kq and we have set 

x = ^K 1 + Q 2 - 2jKQ and Q(t) = qL(t)/(2n). We have also introduced the 



constant J\f 2 defined by 



M 2 



[(f) 
r(l)r(f) 



0.11 for n = 2 , 

0.08 for n = , 

0.05 for n= -3/2 

0.04 for n = -1.8 . 



We have performed the double integrals Ij(k,t) numerically for different values 
of the spectral index n > —2. In the numerical integration we neglect the time 
dependent cutoff in the above integrals and consider only the time dependence given 
by K{t) = K^r 13 , with K* = kLj(2n), r = [(t - t- m )/t* L ] and (3 = 2/(n + 5) for the 
selective decay regime while f3 = 2/3 in the inverse cascade phase. In general we 
find that the antisymmetric contribution J 2 is negative (as it should be |0|), and 



negligible with respect to the symmetric contribution, namely 1\ <C |/ 2 |. 

First we consider a blue magnetic field spectrum which is characterized by n > 
—3/2. The result of the numerical integration for n = 2 , which corresponds to a 
causal magnetic field, is shown in Fig. [□]. The integral I(K) can be approximated 
by the following analytical expression, 



I(K) 



0.034 



i+(sr+(f) 



(B.4) 
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K(t) 

Figure 11: Symmetric part I\{K) (blue, dash-dotted line) and absolute value of the anti- 
symmetric part |/2(-f^)| (green, dotted line) of the anisotropic stresses for a blue magnetic 
field with n = 2. The sum of the two parts I(K) = I\ + AI2 (red, dashed line) and the fit 
(black, solid line) are also shown. 




K(t) 

Figure 12: Symmetric part I\{K) (blue, dash-dotted line) and absolute value of the 
antisymmetric part |/2(iiT)| (green, dotted line) of the anisotropic stresses considering now 
a red magnetic field with n = —1.8. We plot also the sum of the two parts I(K) (red, 
dashed line) and the fit (black, solid line). 

With this approximation the equal time correlator U.g[K,t) can be written as 

U B (K,t) ~ —L 3 (t)p 2 R (t) — . (B.5) 
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For a red magnetic field spectrum with n < —3/2 we find a somewhat different 



oc - , lvlr ^ B , (B.6) 



power spectrum which is shown in Fig. [12| for the case n = —1.8. The integral can 
be approximated by 

i + {K/cy 

where the constants A and B are given by A = 2n + 3 < and B = A + 7 /2 = 
2n + 13/2 and C varies with n and has been chosen to be equal to 1.4 in order to fit 
the numerical result. We recover the asymptotic behavior for the anisotropic stress 
of a red magnetic field found in M , pi 



I(K) ~ K 2n+3 for < 1 , 
I(K) ~ f^ 7/2 for K > 1 . 

This explains also the higher amplitude of I(K) in the red case compared to the one 
obtained for a blue magnetic field: we have an infrared divergence for small values 
of Q. 

Finally for n = —1.8 we find 

U B (K,t) s ^W^ ^g^,, ■ (B.7) 
B.2 GW integrals 

To calculate the integrals (|3.21|) we use Eqs. ( |3.18|) , (]3.19|) for the time evolution of 



the magnetic anisotropic stress and we use the integration variable r 

f - 1- 



T 



t* 



We can then write the GW energy spectrum as 

i2n) 3 K 3 {[I cJ (k)+I CtII (k)} 



dp G w(k,t) _ 2G J\f 2 (p* B ) 2 /o3r/3f[r n , , T n ^ 2 



dlogk 7ra 4 (t) 2n H^Q^ 

+ [I s ,i(k) + I s , n (k)] 2 } , (B.8) 

which yields the following expression for the present density parameter of GWs 

dn GW (k,t ) (fr B ) 2 

— — — ~ 3A/ 2 — Jgw(^*J , 

a log k i 2 rad 

l GW (K,) = K 3 {[I c j(K^ + I cJI (K^] 2 + [I s j(K,) + I S:II (K,)] 2 } . (B.9) 

The four integrals above distinguish the two different phases of the inverse cascade, 
namely the first one where the magnetic energy density is growing linearly up to its 
maximum value p* B (t in < t < £*) and the second one where ps(i) decays as t~ 2 / 3 
(t > Q. 
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For the EW phase transition, n = 2 we have for t- m < t < 



I c j{K*) = / dr 



-f 

Jo 



V0.034 r 2 



x 



12 



6 



3\ 7 / 2 



COS 

• 1 



. A I — • ' 

e 



i s A K *) = / rfr 



o 



sin 



7/2 



X 



-'A I — • - 

2w L e 



(B.10) 



(B.ll) 



while for t >t* the integrals can be written as 



Ic,Il( K *) 



-I 



Is,Il( K *) 



dr- 



V0.034 r 1 / 3 



x 



(T + ^)j 1 + (i^) + ( 



if*T 2 / 3 



7/2 



COS 



1' A I — - - 

vl e 



'1 



V0.034 r 1 / 3 



+ V 1 + W + (^) 



x 



7/2 



sin 



2tt/u I — • - 
2w L e 



(B-12) 



(B.13) 



Secondly we consider a red magnetic field with n = —1.8 < —3/2, that implies 
A = — 3/5, B = 29/10 and C = 1.4. In this case the four integrals read for t m <t<t* 



I cJ {K*) = / dr 



r 9 / 5 (i^/40)- 3 / 10 



x 



cos 



2%K, 



r 1 
2^ + 7 



= / dr 



T 



9 / 5 (if*/40) 



-3/10 







sin 



( T +1^1 +(*£!!) 



X 



29/10 



2tt^ 



r 1 
2^~ + 7 



(B.14) 



(B.15) 
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while for t > t* the integrals can be written as 

r 2 / 15 (^/40)- 3 / 10 



.11 



dr- 



x 



( T + Wl+(** 



.4 



29/10 



COS 



2^ 



(B.16) 



'fin 



-2/15 



(K*/4Q) 



-3/10 



X 



29/10 



sin 



T 



2vt 



(B.17) 



Inserting typical values for the above quantities we perform a numerical integra- 
tion, and we find that the first phase, the 'switching on' of the inverse cascade, is 
completely irrelevant for the final result for most of the spectrum. It does, however 



affect the peak position and the decay law as we discuss in Section |3.4j . The numerical 
solutions of the integrals are shown in Figs. |6] and 

B.3 The fits for the GW spectrum 

In deriving the analytical fits to the numerical GW spectra, Eqs. (|3.30|) and (|3.32| ), 
we have been guided by analytic intuition of the behaviour of the integrals given in 
Eqs. ( B.10|) to flB.lTp above. Here we give some details for the understanding of the 
fits. 

Let us start with the causal case, n = 2. First of all, the main contribution 
to the GW spectrum comes from the integral in Eq. ( p,12|) , i.e. the cosine part in 
( |3.21| ). For very small values of K*, below the characteristic wave number k < l/t& n , 
the cosine does not oscillate: therefore, we expect to inherit directly the slope of 
the anisotropic stress. For the causal case this is flat, consequently we expect a 
Kl behaviour, coming from ( |B.9|) . The constant €\ is fixed by the large wavelength 
limit of Zgw(-K*)j given mainly by the integral of Eq. (|B.12|) evaluated at the upper 
boundary t^ n : 



J GW (/C ->0) ~ 7.73 n(o) 



An 



2/3 



Kl = a Kl 



For higher values of the wave number, the main contribution to the integral comes 
roughly from the first oscillation of the cosine in Eq. (p.l2j ) (note that the integrand 
decays with time). This can be accounted for by integrating only up to the time 
t — l/k, causing a change of slope of the GW spectrum, which now results in 
Igw{K*) oc kV 3, . In the main text this slope is set to which corresponds to 
the best fit result from the numerical evaluation of the integral (see Fig. 13). These 
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analytical considerations are in fact quite crude and lead to slopes which are not 
very precise. The parameter e 2 is determined by the matching at the limiting value 
k = 1/tfln, which is the value of the wave number for which the cosine starts to 
oscillate: 

6 2 ^0.07^Ie^ X 

This behaviour continues until k becomes of the order of Above this value, 

the time dependence of the integrand is no longer r _2//3 but r 1 / 3 (see Eq. (|B.12|) ). 



This results in a further change in the slope of the spectrum, which now becomes 
^gw(^) K^ 3 . In the main text this slope is set to \/TQ, again according to the 
numerical evaluation of the integral. By continuity, the paramter €3 in Eq. ( |3.30| ) is 

/ rp \ 1/3 

In the inflationary case, the main contribution to the GW spectrum comes again 
from the integral in Eq. ( p.l6|) . For very small values of K*, below the characteristic 
wave number k < l/t nn , we expect to inherit the slope of the anisotropic stress. The 
constant €4 is given by the large wavelength limit of Zgw(-K*) : 

/ T \ 4/15 

J GW (#* -> 0) ~ 68 IL(K.) ( ^ ) Kl = e 4 . (B.18) 



fin 



The above formula is valid for n = —1.8. For higher values of the wave number, the 
same argument as in the causal case applies, and we integrate only up to t ~ 1/k: 
the slope in wave number of the GW spectrum now results in Z G w(-^C) ^ iQ 2n+10 ^ 3 . 
By continuity, we obtain (again for n = -1.8) e 5 ~ 619/(2tt) 4 / 15 . This behaviour 
continues until the wave number for which the final time of turbulence tfi n (k), given in 
Eq. ( 3.26Q , becomes smaller than 1/k: this happens for ~ ((27r) 5 /i?^) 1 / 7 , see the 



first line of Eq. (|A.17f) . For higher wave numbers, the upper limit of integration has 
a different A;— behaviour which translates to the slope Xgw(-^*) oc K* ( - 2+6n ^ 5 . By 
continuity, e 6 ~ 619/ Rl 2 ^ 25 ~ 22 for n = — 1.8 and ~ 10 3 for inflation. 
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